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Abstract. This work deals with both instantaneous uniform mixing property and temporal 
standard deviation for continuous-time quantum random walks on circles in order to study their 
fluctuations comparing with discrete-time quantum random walks, and continuous- and discrete- 
time classical random walks. 



1 Introduction 

In this letter we consider mainly temporal fluctuations of continuous-time quantum random 
walks (QRWs). Farlii and Gutmann Ij proposed the model of a continuous-time QRW 
which might be useful for quantum computation. Recently Childs et al. P] showed that an 
exponential speedup was achieved by using continuous-time QRWs. Moore and Russell [3j 
investigated the mixing times and distributions for the continuous-time case on the hypercube 
first. Furthermore, Ahmadi et al. |3] and Adamczak et al. [S] studied mixing properties 
on various graphs, e.g., complete graph, ra-cube, cycle. Gerhardt and Watrous proved 
non-instantaneous uniform mixing on some Cayley graphs of the symmetric groups Sn for 
n > 3. On the other hand, a discrete-time model was introduced by Ambainis et al. [Zj and 
Aharonov et al. , and has been explored by several groups (for examples, j^l UHl UH 1121 

HSHH). 

The present letter treats a continuous-time QRW on C^- which is the cycle with N 
vertices, i.e., = {0, 1, . . . , — 1}. Let A be the N x N adjacency matrix of Cn- The 
continuous-time QRW considered here is given by the following unitary matrix: 

U{t) = e-'*^/2 

Following the paper of Ahmadi et al. |4 , we take —itA/2 instead of itA/2 . The amplitude 
wave function at time t, |\l/Ar(t)), is defined by 

\^N{t)) = u{t)\^Nm 

The [n + l)-th coordinate of |\&Ar(t)) is denoted by |\l/Ar(?T,, t)) which is the amplitude wave 
function at vertex n at time t for n = 0, 1, . . . , — 1. In this letter we take |\l/jv(0)) = 
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[1, 0,0,..., 0]^ as an initial state, where T denotes the transposed operator. The probabihty 
that the particle is at vertex n at time t, P]y{n,t), is given by 

PN{n,t) = {^N{n,t)\^N{n,t)) 

For the details of the definition for continuous-time QRWs are given in jH E] . 

Here we give a connection between a continuum time limit model for discrete-time QRW 
on Z (= the set of integers) given by Romanelli et al. Jl] and continuous-time QRW on 
circles. Ahmadi et al. ^ showed 

1 N-l 

l^iv(n,t)) = -Ee-*(*^°^«^-"«^), (1.1) 
j=o 

for any t > and n = 0, 1, . . . , — 1, where C,j = 271] /N. As Ahmadi et al. |3] pointed out, 
\^N{n,t)) can be also expressed by using the Bessel function, that is, 

\^N{n,t))= i-^fUt) = \ E HM(t), (1.2) 

k=n imodN) k=±n (modN) 

where Jv{z) is the Bessel function (see |il5j). In fact, combining the generating function of 
the Bessel function: 



exp 



and Jy{—z) = {—iyjiy{z) gives 



-t 













By using Eq. (jl.lj) . the amplitude wave function is then calculated as 

1 Af-l 

= T7 E E(-0^e^^"+'')«^ J.(t). (1.3) 
j=o i/ez 

In a similar way, we have 

1 N-l 

\^N{n,t)) = - E E(-0^e^^"-''^^^ J.(t). (1.4) 
From Eqs. ()1.3|) and p.4|) . we obtain Eq. p.2|) . Moreover, Eq. ()1.2|) gives 

= E JkN+nitf + E '^OS (^(i - JjN+n{t)JkN+n{t). (1.5) 

The following formula is well known (see page 213 in Ref. |15j): 

E Mtr = 1, (1.6) 

k£Z 
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for t > 0. That is, {Jfc(t)^ : A; G Z} is a probability distribution on Z for any time t. Eq. 
()1.6p can be rewritten as 



iV-l 

E E Wn(t)' = 1, (1.7) 
n=0 fcGZ 



for t > 0. Eqs. dUS)) and (III7|) give 

/vr \ 

^ ^ cos f -(j - /i;)A^j JjN+n(t)JkN+n(t) = 0. 
n=0 j^k:j,k&Z ^ 

Very recently, Romanelli et al. |jl4j studied a continuum time limit for a discrete-time 
QRW on Z and obtained the position probability distribution. When the initial condition 
is given by a;(0) = 5i,0)^K0) = in their notation for the Hadamard walk, the distribution 
becomes the following in our notation: 

P(n,t) = Jn{tl^f. 

More generally, we consider a discrete-time QRW whose coin flip transformation is given by 
the following unitary matrix: 



a b 
c d 



Note that a = b = c = —d = l/\/2 case is equivalent to the Hadamard walk. In a similar 
fashion, we have 

P{n,t) = Jniatf. (1.8) 

In fact, Eq. (jl.fjj) guarantees J2nez Pip-^ ^) = 1 ^"^J t > 0. From the probability theoretical 
point of view, it is interesting that a continuum time limit model for a discrete-time QRW 
on Z gives a simple model having a squared Bessel distribution in the above meaning. For 
large A^, Eq. ()1.5|1 is similar to Eq. ()1.8|1 except for the second term of the right-hand side 
ofEq. |T3|1 . 

To investigate fluctuations for continuous-time QRWs, we consider both instantaneous 
uniform mixing property (see below) and temporal standard deviation (see Section 4). The 
former corresponds to a spatial fluctuation and the latter corresponds a temporal one. 

A random walk starting from on C^v has the instantaneous uniform mixing property 
(lUMP) if there exists t > such that P7v(n, t) = 1/A^ for any n = 0, 1, . . . , - 1. To restate 
this definition, we introduce the total variation distance which is the notion of distance 
between probability distributions defined as 

\\P-Q\\ = m^x\PiA)-Q{A)\ = l Y \P{n)-Qin)\ (1.9) 

^^^^ ^ neCjv 

By using this, the lUMP can be rewritten as follows: if there exists t > such that | \PN{-,t) — 
71n\ \ = 0, where hn is the uniform distribution on Cn- As we will mention in the next section. 
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continuous-time QRWs on Cat for = 3,4 have the lUMP. On the other hand, Ahmadi et 
ah |lj conjectured that continuous-time QRWs on C^v for A^ > 5 do not have the lUMP. It 
will be stated in the last section that the continuous-time classical random walk on Cn does 
not have the lUMP for any A^ > 3. In our previous paper JH], we consider the temporal 
standard deviation for discrete-time QRWs on Cn, in particular, when A^ is odd. Similarly 
we study it here for continuous-time case and clarify some differences between both cases. 
The definition of the temporal standard deviation will be given in Section 4. 

The structure of the present letter is as follows. Section 2 is devoted to a review on 
QRWs on circles. In Section 3, we consider probabihty distributions and time-averaged dis- 
tributions for QRWs on Cn- Section 4 gives our new results on temporal standard deviation 
on continuous-time quantum case comparing with discrete-time quantum case and classical 
case. In the last section, we discuss a future problem on convergence. 



2 Classical Case 

First we review the classical case on Cat, (see Schinazi ^7] and references in his book, for 
example). In the continuous-time classical random walk, 

1 7V-1 
j=0 

for any t > and n = 0,1, . . . , N — 1. In this walk, the distribution for the random time 
between two jumps is the exponential distribution. As for the discrete-time classical random 
walk, 

PN{n,t) = j-jYl cos(^jn)(cos^j)*, 

j=0 

for any t = 0, 1, ... , and n = 0,1, . . . , N — 1. Here a time-averaged distribution in the 
continuous-time case is given by 

Pjv(n) = lim ^ r PN{n,t)dt, (2.10) 

if the right-hand side of Eq. ()2.10|1 exists. As for the discrete-time case, 

P^(n)= hm ^J2Mn,t), (2.11) 

if the right-hand side of Eq. ()2.1H) exists. Concerning the continuous-time classical case, 
it is easily shown that P/v(n, t) 1/N {t oo) for any n = 0, 1, . . . , A^ — 1. Then we 
have immediately -P/v(n) = 1/A^ for any n = 0, 1,...,A^ — 1. On the other hand, in the 
discrete-time classical case, if A^ is odd, then PN{n,t) 1/N as t — oo, and, if A^ is even, 
then it does not converge in the limit of t ^ oo. However, an important property is that P/v 
is uniform distribution in both continuous- and discrete-time classical cases. 
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3 Quantum Case 

Now we consider continuous-time QRWs. As we will see later, in contrast to classical 
case, P/v does not become uniform distribution for any > 3 in quantum case. From Eq. 
(jl.lj) . we have 

PN{n, ^) = ^ + j;^2 ' (3-12) 

where 



RN{n,t) = cos{t(cos^j - cos^fc) - n{^j - C.k)}, (3.13) 

o<i<fc<Af-i 

for any t > and n = 0, 1, . . . , — 1. It is easily checked that 

RN{n,t) = Rn{N -n,t), 

for any t > and n = 1, . . . , N, where N = [{N — 1)/2] and [x] is the smallest integer greater 
than X. For example, when A^ = 3, 

i?3(0,t) = 2cos(3t/2) + 1, 

i?3(l,t) = i?3(2,t) = -cos(3t/2)-l/2, 

when A^ = 4, 

i?4(0, t) = cos(2t) + 4 cos(t) + 1, 
R^il, t) = i?4(3, t) = - cos(2t) - 1, 
i?4(2, t) = cos(2t) - 4 cos(t) + 1. 



So if A^ = 3, then R3{n,t) = for any t = ±Att/9 + Amn/3 (m G Z) and n = 0, 1,2, and 
if A^ = 4, then R4^{n,t) = for any t = 7r/2 + mvr (m G Z) and = 0, 1,2,3. Therefore 
continuous-time QRWs on C3 and C4 have lUMP. This result was given in Ref. |3]. On the 
other hand, when A^ = 6, 

i?6(0,t) = cos(2t) +4cos(3t/2) +4cos(t) +4cos(t/2) + 2, 

R6{l,t) = i?6(5,t) = -cos(2t) -2cos(3t/2) -cos(t) + 2cos(t/2) - 1, 

i?6(2,t) = i?6(4,t) = cos(2t) - 2cos(3t/2) + cos(t) - 2cos(t/2) - 1, 

i?g(3,t) = -cos(2t) +4cos(3t/2) -4cos(t) -4cos(t/2) + 2. 

Direct computation implies that a continuous-time QRW on Cg does not have lUMP. 

Now we consider the time-averaged distribution for continuous-time case. By using Eq. 
dSini), we have 

P„(«)4 + ^. (3.14) 
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for any n = 0, 1, . . . , — 1, where 



if = odd, ^2n 7^ (mod 2it), 
RN{n) = { -I if N = even, ^2n 7^ (mod 2n), 

if 6^ = (mod27r). 




This result gives 



As for discrete-time quantum case given by the Hadamard transformation, Aharonov et 
al. P and Bednarska et al. [12] showed that Pwin) = 1/N if N is odd or A^ = 4, and 
P/v(n) ^ 1/A^, otherwise. So in contrast to classical continuous- and discrete-time and 
quantum discrete-time cases, P/v in quantum continuous-time case is not uniform distribution 
on Cn for any A^ > 3. In fact, when A^ = odd (i.e., N = {N - l)/2). 



N-l 



p - A^-ll 1 1 Ix 

when N = even (i.e., N = {N - 2)/2), 

{N-2)/2 {N-2)/2 



_ 1 A^-2 1 2 1 2 1 A^-2 1 2 1 



^A^ ' N N^'"''N N^' N ' N m'"''N 

For examples, 

5 2 2\ ^ /3 1 3 1 



9'9'9/' V8'8'8'8 



4 Temporal Standard Deviation 

We consider the following temporal standard deviation aN{n) in the continuous-time case 
as in the discrete-time case: 



a^{n)=\im^^^ [P^M - PN{n)f dt, (4.16) 
if the right-hand side of Eq. ()4.16|) exists. In the discrete-time case. 



1™, ^ f E {PNin, t) - PN{n))\ (4.17) 



if the right-hand side of Eq. ()4.17|) exists. In both continuous- and discrete-time classical 

cases, 

(Tiv(n) = 0, 



6 



for N > 3 and n = 0, 1,...,A^ — 1. The reason is as follows. In the case of classical random 
walk starting from a site for = odd (i.e., aperiodic case), a coupling method implies that 
there exist a e (0, 1) and C > (are independent of n and t) such that 
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\PN{n,t) - PN{n)\ < Ca' 
for any n and t (see page 63 of Schinazi jTTj , for example) . Therefore we obtain 

iE(P.(n,t)-P.(n)) <^^. 

The above inequality implies that o"Ar(n) = (ri = 0, 1, . . . , — 1). As for A^ = even (i.e., 
periodic) case, we have the same conclusion (TN{n) = for any n by using a little modified 
argument. The same conclusion can be extended to continuous-time classical case in a 
standard fashion. 

For discrete-time quantum case, our previous paper [THj implies 

al(n) = ^ [2 {Slin) + SUn)} + US^o + lO^o^i + 35^ - S^in)] - (4.18) 



for any n = 0, 1, . . . , A^ — 1, where A^(> 3) is odd and 

^-^ 1 ^-^ cose- 

s = 'V^ s = \^ 

3 + cos 9j ' 3 + cos 9j ' 

^^ cos((n-l)gj) + cos (ne,) 
^^^""^ - j^, 3T^, ' 

_ cos {{n - 1)9 j) - cos {nOj) 
^^""^ ~ 3 + cos^, ' 

^ , , 7 + cos(2^j) + 8 cos cos2 {{n - 1/2) 6/^ 

^^^"^ - h (3 + cos^,)^ 

with 6j = = 4:7rj/N. For example. 



2V46 2 
as{0) = as{l) = ^, ^3(2) = -. (4.19) 



Furthermore we ^HI obtained 



. 13 - 8v^ / 1 , , , 

M0)= ^ +0 - (4.20) 



A^ VAT 

as A^ — >• oo. The above result implies that the temporal fluctuation o"iv(0) in discrete-time 
case decays in the form 1/N as N increases. 
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Now we consider continuous-time quantum case. For example, direct computation gives 

2^2 



0^3(0) 
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(74(0) = (T4(2) 
fX6(0) = (T6(3) 



^3(1) = ^3(2) 

^34 



9 ' 



16 



0-4(1) = 0-4(3) 



V2 
16' 



2V2 



From now on we focus on general N = odd case to obtain results corresponding to 
discrete-time case (i.e., Eqs. ()4.18j) and ()4.2()j) ). The definition of a^in) implies 



4 
4 

7 X 

iV4 I T- 



(^Nin) = J^^^^_I^fJ^ {RN{n,t) - RN{n)f dt 
lim ^ r RN{n,tfdt-RN{nf 

T^oo T JO 



For = odd case, we have 

i?Ar(n) = 

From Eq. ()3.13p . we have 
1 /"^ 

lim - / R%{n,t)dt 

T-»oo 1 Jo 



{N-l)/2 n = 0, 
-1/2 



n= l,2,...,Ar-l. 



lim — / 

T^oo T Jo 



X 



J2 cos{t(cos^j, - cos^fcj - n{^j^ - 61)} 

0<ji<fci<Af-l 

J2 cos{t(cos^j2 - cos^fcj - n{^j.^ - ^^J} 

0<j2<k2<N-l 



dt 



= — r dt y y 

^00 z J .70 o<ji<fci<Af-l 0<j2<fc2<A'-l 

[cos{t((cos^j, - cos^fej - (cos^j2 - cos^fcj) - n{{^j^ - ^^J - (^^^ - ^^J)} 

+ COs{t((cOS^j, - COS^fcJ + (cOS^j2 - COS^fcJ) - n((Cj, - ■^fci) + {^32 - -^fca))}] 



Remark that each (A^^ — A^)^ term becomes zero except the following two cases: 

COS^j^ — COS,^fc^ — COS^j2 + COS^fc2 = 0, 
COS^jj — COS^fc^ + COS^j2 ~ COS^fcj = 0. 
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Therefore some complicated computations imply 



\\m RUO,t)dt = |(iV-l)2 + l(iV-l)(5iV-13) 

T^oo J Jo o O 



l(iV-l)(5iV-9) 



and 

I „o . . , 1 



lim- / Rl(n,t)dt = -l2N^-AN + 6 — 

T^ooT Jo ' 8 V cos2^ 



+\ ( -6iV + 12 + 



cos2^^ 



for any n = 1, 2, . . . , — 1. Then we have the following main result: when A^(> 3) is odd, 



^^(0) = ^^ Jf^ ' ^^(^) = ^^ (n=l,...,iV-l). (4.21) 

In particular, crAr(O) > a^in) for any n = 1, . . . , — 1. We should remark that there is a 
difference between continuous-time case and discrete-time one, since, for example, Eq. ()4.19|1 
gives (73(0) = (73(1) > (7^(2) in discrete-time case. Furthermore Eq. ()4.21|) implies 



a^(0) = - + o(^-j, a^(n) = - + o(^-j (n = 1, . . . , iV - 1) (4.22) 

as 00. The above result imphes that the temporal fluctuation crAr(O) in continuous- 

time case decays in the form 1/N as N increases as in the discrete-time case. However, from 



13 — 8-\/2 < 2 (see Eqs. fl4.2()|l and ()4.22j) ). it is obtained that a temporal fluctuation of 
continuous-time case at site is greater than that of discrete-time case for large N . 



5 Future Problem 

Finally we consider a future problem on the total variation convergence and weak conver- 
gence. We should remark that the total variation convergence is stronger than weak con- 
vergence. Parity problem vanishes in continuous-time case, so we focus on continuous-time 
case even in the classical case. Now we introduce 

C?Ar(t) = ||PAr(-,t) - TTAfll = i ^ |PAr(n, t) - TTat (n) | , 

where tin is the uniform distribution on Cn- In the classical continuous-time case, dMit) 
decreases as t increases, and 

lim dj^it) = 0. 



9 



So the continuous-time classical random walk on Cat does not have the JUMP for any > 3. 
Moreover, it is known that 

lim dN{tN^) = dooit) = 11' \ftix) - l\dx, 

N-*oo Z Jo 

where ft is the density of mod 1 and Z has the standard normal distribution (see 
Chapters 2 and 5 in pfH]). This fact corresponds to weak convergence of classical random 
walk to Brownian motion on a torus T, i.e., R mod 1. In the continuous-time quantum 
case, an interesting future problem is to obtain a limit theorem corresponding to the above 
classical case. 
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